Evidence for magnetic clusters in BaCoOs 
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Magnetic properties of the transition metal oxide BaCoOs are analyzed on the basis of the ex- 
perimental and theoretical literature available via ab inito calculations. These can be explained by 
assuming the material to be formed by noninteracting ferromagnetic clusters of about 1.2 nm in 
diameter separated by about 3 diameters. Above about 50 K, the so-called blocking temperature, 
superparamagnetic behavior of the magnetic clusters occurs and, above 250 K, paramagnetism sets 
in. 

PACS numbers: 71.15.Ap, 71.15.Mb, 71.20.Be, 75.20.-g 
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Transition metal oxides are among the most stud- 
ied materials in recent literature. The variety of their 
physical properties has made them a topic of inten- 
sive research during the last few years, the phenom- 
ena of colossal magnetoresistancei and high temper- 
ature superconductivity 2 being the most remarkable. 
Among them, Co oxides have drawn much attention, 
since they posses both types of phenomenal. The mag- 
netic properties of such oxides are a matter of ongo- 
ing work. The appearance of properties such as orbital 
ordering^, charge ordering^, spin glass behavior or phase 
segregationiSiSiifliii makes them very rich in phenomena, 
whose understanding is a challenge for the scientific com- 
munity. 

BaCoOs is a transition metal oxide that recently has 
been a matter of interest of both, experimental and the- 
oretical WO rklSii2ildti5iiStil. Its structure^ can be de- 
scribed as a 2H-hexagonal pseudoperovskite formed by 
chains of distorted, tilted, and face-sharing Co06 oc- 
tahedra. The plane perpendicular to the chains con- 
sists of a hexagonal array of Co atoms, which would 
lead to frustration for any in-plane collinear antiferro- 
magnetism. The only possible long-range antiferromag- 
netic state (for collinear moments) could couple ferro- 
magnetic planes antiferromagnetically (so-called A-type 
structure The magnetic properties of BaCoOs have 
still not been resolved in literature. It is certain that the 
Co 4+ ions are in a low-spin state (S=l/2, t| ff e°), since 
an ion with such a high valency produces a large crys- 
tal field. This has been shown both experimentally 12 
and theoreticalljii 4 * 1 ^* 1 ^. Moreover, the experimentally 
found semiconducting behavior— has been predicted 
theoreticalljiii, using the LDA+U approach, to be a con- 
sequence of an orbital ordering phenomenon, which is not 
yet confirmed experimentally. The actual magnetic con- 
figuration of the system is still uncertain. So far no neu- 
tron diffraction measurements are available in literature 
and the experiments that have measured the suscepti- 
bility of the systemic are not conclusive. Theoretical 
APW+lo LDA+U calculations^ predict a ferromagnetic 
state as the magnetic ground state. The most stable 
antiferromagnetic state is higher in energy but this en- 



ergy difference is extremely small and depends on the 
value of U. Therefore any magnetic coupling that might 
occur will have a very small stabilization energy. This 
feature is typical for systems showing a cusp in their 
susceptibility^, as for BaCoOs— 

In this paper, we try to shed some light in the explana- 
tion of the magnetic properties of BaCoOs and resolve an 
open question in literature. As a start we will analyze the 
experimental data available and propose a model based 
on the existence of noninteracting ferromagnetic clusters 
in the system, which can explain the observations. Then, 
using ab initio calculations, we will predict the size and 
density of the clusters, establishing a picture which is 
fully consistent with both experimental and theoretical 
data. 

When we take a close look at Ref. 0, Fig. 10, we 
observe susceptibility curves that are typical for a "fine- 
particle" system^ with a blocking temperature of around 
50 K. This is approximately, where the field cooling (FC) 
magnetization curve separates from the zero-field cooling 
(ZFC) one for an applied field of 1 kOe. These FC and 
ZFC curves indicate that BaCoOs is formed by regions 
or clusters, whose magnetic moments are ferromagneti- 
cally ordered, but are dispersed into a non-ferromagnetic 
matrix, a situation that can be described by Wohlfarth's 
superparamagnetic modelSi. The shape of the FC curve 
below the blocking temperature indicates that the clus- 
ters are widely separated from each other and they do 
not interact strongly. However, would they interact, this 
curve would be flat instead of growing exponentially as 
is observed— 

The blocking temperature depends on the size of the 
clusters, the anisotropy constant of the material, the ap- 
plied field and on the measuring time (on the apparatus 
utilized for the experiment). By analyzing the data from 
Ref. [L2J, one can see that the blocking temperature de- 
creases when the applied field increases. This is the typ- 
ical behavior of fine-particle systems, even though the 
contrary might occur if the clusters grow with the ap- 
plied field 22 . Based on the theory developed by Neel for 
superparamagnetic particle o 23 i 24 , we can predict that the 
highest blocking temperature (which is reached when the 
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applied field is much smaller than the so-called anisotropy 
field) is basically obtained at IkOe. According to this, 
the mean blocking temperature varies as: 
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where Tb is the blocking temperature and Hk is the 
mean anisotropy field. From the data in Ref. lL2l 
the blocking temperature is approximately 10 K for 
H=10kOe and 50 K for H=lkOe. From this, the maxi- 
mum blocking temperature that can be reached at a very 
small field is about 53 K, a limit that is basically obtained 
at a field of IkOe. This is caused by the large anisotropy 
of the system, which is due to the quasi-one-dimensional 
structure formed by widely separated chains of Co atoms. 
Hence, the value we have chosen to explain the magnetic 
properties of the system is Tb— 50 K, the expected value 
for low fields. The description of the system, which will 
be given below, corresponds to these low applied fields. 

From the value of the blocking temperature we can 
estimate the mean size of the clusters. For doing so, 
we have calculated the anisotropy constant of the sys- 
tem from ab initio methods. The calculations were per- 
formed using the WIEN2k2^2& software, a package that 
uses a full-potential, all electron APW+lo2£ method that 
allows to carry out total energy calculations with one 
of the most precise methods available. For this mod- 
erately correlated transition metal oxide, we used the 
LDA+U approach including self-interaction corrections^ 
in the so-called "fully- localized limit"— with U= 5eV and 
J= 0.5eV, values discussed in Ref. [l7. This method 
has proven reliable for transition metal oxides 3 ^, since 
it improves over GGA (generalized gradient approxima- 
tion) or LDA (local density approximation) in the study 
of correlated electrons by means of introducing the on- 
site Coulomb repulsion U. The non-orbital dependent 
part of the exchange-correlation potential was calculated 
using the GGA in the PBE (Perdew-Burke-Ernzerhof) 
schemed. Local orbitals (Co 3s, 3p, 3d, Ba 4d, 5s, 5p, 
O 2s and 2p) were added to obtain a better flexibility 
of the basis set and to improve the description of the 
semicore states. R m tK maa; =7 and 500 k-points in the 
irreducible Brillouin zone (1200 k-points in the whole 
Brillouin zone) were taken. For calculating the magnetic 
anisotropy, a fully relativistic calculation was carried out, 
where spin-orbit effects were included using a second vari- 
ational scheme 3 ^. 

For calculating the energy contribution of the magnetic 
anisotropy we studied the system in different configu- 
rations with the magnetic moments along the different 
crystallographic directions: (100), (110) and (111), and 
compared their total energies. The system was studied 
in a ferromagnetic low-spin state with an alternating or- 
bital order along the c axis as described in Ref. IT7L which 
in this paper was shown to be the ground state. From 
our total energy calculations, the value of the anisotropy 
constant is 1.2 mRy per unit cell, which contains 2 Co 



atoms (about 2xl0 8 erg/cm 3 ). This value is larger than 
the anisotropy energy of a system of fine Co particles 
(~ 7xl0 6 erg/cm 3 ). We are not considering the shape 
anisotropy due to the formation of clusters, but we can 
assume this contribution to be a secondary effect in this 
highly anisotropic compound. The lowest energy state 
has the moments lying in the hexagonal plane (along the 
b axis). We assume the system to have cubic anisotropy 
(since the Co ions are in an octahedral environment) and 
found the value of the anisotropy constant as K = ^f-, 
where Eg is the energy difference between the easy (b) 
and hard (c) axis. The procedure for estimating the vol- 
ume of each cluster is via the following formula for su- 
perparamagnetic particles^ 3 -: 



KV = 2bk B T B 



(2) 



where Tb is the blocking temperature, V is the mean 
volume of the cluster and kg is the Boltzmann constant. 
The factor 25 comes from the measuring time, which was 
assumed with a time constant r= 10 2 s for the experi- 
ments of magnetic susceptibility under FC and ZFC con- 
ditions in Ref. ^3- Using equation [2 with the anisotropy 
constant calculated ab initio and the blocking tempera- 
ture from the experimental data in Ref. [12J (Tb — 50 
K), we estimate that each cluster has a typical diameter 
of approximately 1.2 nm, containing about 14 Co atoms. 
The value of the anisotropy constant was calculated at 
T= 0, hence the volume of the cluster must be considered 
as a lowest limit. There exists a distribution of particle 
sizes. This is evident by the fact that the maximum of 
the ZFC curve is displaced from the point where ZFC and 
FC curves separate from one another. In our description, 
we will stick to the average values. 

From the blocking temperature (at about 50 K) to ap- 
proximately 250 K, the system shows a superparamag- 
netic behavior, entering the normal paramagnetic regime 
at about 250 K. This superparamagnetism can be due to 
the presence of magnetic clusters, which can be identi- 
fied as an assembly of single domain particles with a total 
magnetic moment of about 14 iir per particle (from the 
ab initio calculations in Ref. Il7l the total magnetic mo- 
ment per Co atom is 1.00 /is), where the magnetization 
of the sample follows the Langevin formula: 



M = N[iL{y) 



(3) 



where N is the density of clusters per unit volume, fi is 
the magnetic moment of each cluster, y — -^M^ and L(y) 
is the Langevin function. At small y, the magnetization 
varies linearly with H/T and the Langevin function ap- 
proximates to a Curie law: 



M ~ N 



3k B T 



(4) 



Figure ^ shows the data from Ref. as M vs. 
H/T. This gives a straight line, that corresponds to 
the Langevin formula (superparamagnetism) at "normal" 
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FIG. 1: Total magnetization vs. H/T for four values of the 
temperature (data from Yamaura et al*—). In the inset, the 
curves at T= 150 K and T= 200 K are highlighted, showing an 
excellent H/T superposition. However, for T< 50 K, straight 
lines have a smaller slope due to the anisotropy effects at 
temperatures comparable to the blocking temperature. 



temperatures (50<T(K)<200) and "moderate" fields 
(0<H(Oc)<5xl0 4 ). But, the same type of plot has a dif- 
ferent slope for T= 50 K and T= 20 K, because they are 
below the blocking temperature and the clusters are said 
to be "blocked" (the effects due to the blocking are appre- 
ciable in a region around the blocking temperature, su- 
perparamagnetism is observed where ZFC and FC curves 
match each other, above some 70 K, from the curves in 
Ref. Il2|) . The slope increases with T and approaches a 
maximum when T >> Tb, being the difference remark- 
able in the vicinity of T#, where thermal equilibrium 
disappears. 

From these curves we have estimated the density of 
clusters that exists in BaCoOa, being the result that we 
have one 1.2 nm-diameter cluster in a spherical volume 
with a diameter of about 3.0 nm. Hence, the average 
distance between two clusters is about 3 diameters, big 
enough to assume that they do not strongly interact with 
each other. As we mentioned above, this fact coincides 
with the shape of the curves in Ref. 12 that ressemble 
noninteracting superparamagnetic particles 20 . 

The normal paramagnetic state appears above some 
250 K, when thermal energy is enough to break the mag- 
netic order of the crystal. Below the blocking tempera- 



ture, Tb — 50 K, magnetic moments within the cluster 
lie along the easy axis and for this system, formed by 
ferromagnetic clusters, hysteresis loops are expected. 

One of the possible mechanism for the formation of 
these clusters is Nagaev's theory^ 4 *^, according to which 
a conduction electron might become trapped and po- 
larizes ferromagnetically the antiferromagnetic medium 
giving rise to the occurrence of ferromagnetic particles 
dispersed into a non-ferromagnetic matrix, producing a 
phase separation phenomenon. 

In this paper, we have tried to explain the magnetic 
properties of the transition metal oxide BaCoC>3 utiliz- 
ing ab initio full potential APW+lo LDA+U calculations 
to interpret the experimental data available. The theo- 
retical prediction of ferromagnetism as the ground state 
does not coincide with the experimental data of mag- 
netic susceptibility, which is not that of a long-range fer- 
romagnet. The explanation we propose is based on the 
formation of noninteracting ferromagnetic clusters in the 
system. Using ab initio calculations, we have estimated 
the distribution of these clusters. They have a typical size 
of about 1.2 nm in diameter, and are separated by about 
3 diameters, which explains the evidences that they do 
not interact strongly. With this model, the susceptibility 
curves can be explained, corresponding to a superpara- 
magnetic behavior in the temperature region from the 
blocking temperature (Tg ~ 50 K) up to 250 K. Also, the 
change in the slope of the magnetization vs. H/T curves 
below 50 K is related to the blocking of the superpara- 
magnetic state. Above 250 K, normal paramagnetism is 
found. 

From our conclusions, BaCoOa is expected to be an ex- 
ample of a system with a similar behavior to that formed 
by "ultra-fine" noninteracting magnetic particles and fur- 
ther experimental work is needed to provide more insight 
into the characterization of the material. SANS (small 
angle neutron scattering) measurements and in addition 
a better static and dynamic study of the variation of 
the magnetization with temperature and applied field is 
strongly encouraged. 
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